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Abstract 

A mean field SDW analysis of pseudogap in the underdoped cuprates is 
proposed on the basis of the t — t' — U Hubbard model. The prediction 
of our theory is consistent with the experiment quite well within the 
uncertainty of the present experimental measurement. Therefore we 
argue that the pseudogap phenomenon in the underdoped cuprates can 
be well explained within the mean field approximation. 
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1. Introduction 



A large body of experimental investigations have indicated that underdoped 
high-temperature superconductors exhibit intriguing properties at temperatures above 
the superconducting transition temperature T c . Most notably, the underdoped cuprates 
exhibit a pseudogap behavior below a characteristic temperature T* which can be 
well above the superconducting transition temperature T c . The so-called "pseudogap" 
means a partial gap. "An example of such a partial gap would be a situation where, 
within the band theory approximation, some regions of the Fermi surface become 
gapped while other parts retain their conducting prope rties and with increased 
doping the gapped portion diminishes and the materials become more metallic" (quoted 
from Ref. [p]]). What is the origin behind it? A number of scenarios like pair 
formation well above T c @-fU, spin-charge separation spin-density wave (SDW) 
or antiferromagnetic fluctuations [0,|J have been proposed as possible origins of these 
pseudogap phenomena. However, no consensus has been reached so far, which one is 
correct in these microscopic theories. It should be noted that these theories of the 
pseudogap are all beyond the mean-field approximation. 

In present paper, we propose a mean field SDW analysis for understanding the 
pseudogap phenomena in the underdoped cuprates. Our aim is to examine to what 
extend can we interpret this phenomenon within mean field theory. 

2. Electronic band structure 

Soon after the discovery of the cuprate superconductors, the electronic band 
structure of the cuprates has been calculated by the local density approximation (LDA) 
band calculation [PHUI . The result of the electronic band structure of the cuprates of 



the LDA band calculation is consistent with the later angle resolved photoemission 
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experiment [0-4T6|. The electronic band structure of the cuprates can be well fitted 



by a tight-binding model, which is written as 

£k = — 2t(cos k x + cos k y ) — At' cos k x cos k y . (1) 

Where t is nearest-neighbor, if is next-to-nearest-neighbor. In this paper we consider 
t > and if < only. Energy contour lines for the electronic band structure ([I]) are 
shown in Fig. |l|. There are two different saddle points locate at the M points [(±7r, 0) 
and (0, ±7r)] of the Brillouin zone. The energy contour line with energy e s = At' pass 
through the saddle points. 

For convenience, we choice M as the new origin of the k-space and take the energy 
e s = At' at the saddle point M as zero. Then the dispersion ([!]) is reexpressed in the 
form 

£k = - At' 

= —2t(— cos k x + cos k y ) + 4t'(cos k x cos k y — 1). (2) 

If without specific statement, we keep this usage later. 

We replot in the period Brillouin zone the energy contour lines passing through the 
saddle points. As shown in Fig. ^ there are two different regions: I + I' and II. In the 
region II, < 0, in the regions I + I', £k > 0. The area of the region I + I' is larger 
than that of the region II. When the region II is shifted by the vector Q = (71*, n), it 
coincides with the region I. The region I' is called as the necklace region, which has 
following features. Firstly, when k locates in a bubble, k + Q will locate in another 
one, both and £k+Q are larger then zero. On the other hand, in the regions outside 
the necklace region I', both sign of the e k and £k+Q are always opposite. For example, 
when k locates in I, e k > 0, then k + Q will locate in II, £k+Q < 0. Secondly, in the 
overdoping regime, the Fermi surface entirely lies outside the necklace region (as shown 
in Fig. |^). But for the underdoping case, only part of the Fermi surface lies outside the 



necklace region, and further, with decreased doping the portion outside the necklace 
region increases (as shown in Fig. |^). 

It is interest to note the fact that when if — 0, the necklace region and said 
peculiarity of the band structure of the cuprates disappears. 

3. Mean-field theory 

The starting point of our calculation is the Hubbard model. In the momentum 
representation, the t — if — U Hubbard model can be written as |TJJ 

H = H( £ k - H)al a a ka - ^- Yl 4+qe^ a kc r 4'cX a k'+q CT '- (3) 

k(T /iV q kcrk'cr' 

Here a term \NU has been omitted. U is the local Coulomb repulsion. a\ S(r (al La .) is 
the annihilation (creation) operator for the electron with momentum k and spin a. 
fi is the chemical potential, e k is given by Eq. (0). All the momentum summations 
extend over the Brillouin zone. Considering commensurate SDW state and using the 
mean-field approximation, the Hamiltonian reduces 

H = J2'( £ * ~ v) a L a ^ + H'^k+Q - /i)4 +Q(T a k +Qcr - A£) {al +Qa aa kc + h.c). (4) 

kcr kcr kcr 

Here means that the sum extends over the magnetic Brillouin zone (shown in Fig. [| 
by the thick square). The term ^jA 2 has been omitted. The order parameter A is 
given by 

2,U i 

A = < 4+Qa^«kcr > • (5) 

kcr 

By the following canonical transformation 

«kcr = WkGtkcr - ^kO"Ctk+Qcr , 

7ka = v k aa ka + u k aa k+Qa , (6) 
the Hamiltonian (ffl) is diagonalised as 
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H = H'(^i(k)aL«k CT + £ 2 (k)7L7k<x), (7) 



in which, 



ei (k) = £k + " k+Q - M + J( £k ~ £k+Q )^ + A 2 , 



ea(k ) = £k + £k+Q — fi — J( £k " £k+Q ) 2 + A 2 , (9) 



-E'4r(tanh(^) - tanh(^l )) (10) 



and 



^(k) = ^/( £k 2 " k+Q ) 2 + A 2 . (11) 

Here £i(k) and £ 2 (k) are energy dispersions of the quasiparticles. For the hole doping 
system, the Fermi surface lies inside the lower band (e 2 (k)). The pseudogap is given 
by 

Ap 5 (0) = h(k)| 



= /i — 4t'(cos k x cos k y — 1) + J 4t 2 (cos k x — cos k y ) 2 + A 2 . (12) 

In Fig. [5] we plot the part of the magnetic Brillouin zone of the Fig. |4j. The light 
curve represents the Fermi surface. k x - and fc y -axis are parallel with Mr and MX, 
respectively. In Eq. (0), k = (k x ,k y ) is the wave vector of the Fermi surface, i.e. 
£k — jj, = 0. <p = arctan(/ca;/ k y ) is polar angle of the wave vector k. For convenience, we 
take <$' = arctan(/ca; / (tc — k y )) as variable instead of the <ft m the following calculations. 

4. Results 

In this section, we analyse the angular dependence of the pseudogap Aps((j>') along 
the Fermi surface. Ouing to the symmetry of the energy spectrum £ 2 (k), our analysis 
can be limited only in the interval < <p' < 45°. 



By solving Eqs. @, ([TOD, (|TTJ) and (|12"D numerically, we compute Aps(4>') 
at T = K in the underdoping regime (/x > 0). In the computation, we choose 
t = 430 meV, t'/t = —0.18, U/t = 0.8 and the hole doping concentration x=0.13. The 
results are plotted as Aps(<fi') versus 4>' curve in Fig. ||. It shows that there is strong 
angular dependence of the Apg(<f>'), as one moves along the Fermi surface from 0' = 
(i.e. near the saddle point M, or at the hot spot) to <p' = 45° (i.e. cold spot). At 
first, we see the maximum pseudogap at 0' = 0. As we go into the necklace region, 
the pseudogap drops quickly and, at approximately 18°, drops down to 2 meV. And 
then, the pseudogap decreases monotonously to Aps(45°). Experimental measurement 
reveals that only a portion of the Fermi surface near the saddle point M becomes 
gapped while in other parts, the pseudogap is equal to zero fiyi8fl. However, the 
error-bar of the pseudogap data is rather larger[f. It is impossible to say certainly 
that along the part of the Fermi surface near the cold point, the pseudogap is real 
zero or only a small quantity. Keeping this fact in mind, we conclude that the general 
structure of the pseudogap along the Fermi surface, shown in Fig. [], captures the main 
feature of experiment [p~8|fp. 

The dependence of Aps(0) on the hole doping concentration is shown in Fig. [7]. It 
shows that Aps(0) increase with the decrease of hole doping. In Fig. [|, we plot A(f)' 
versus the hole doping concentration curve. Here, A<fi' is the interval <$' (measured 
from (J)' = 45°), defined by the requirement that the value Aps(4>') is less than a proper 
chosen value (say, 2 meV in Fig. |6| and §). It can be seen from Fig. |]that the length of 
the Fermi arc, along which the pseudogap less than 2 meV, increase with the increase of 
doping. It implies that as doping increase, the portion of the Fermi surface destroyed 



1 ^See, for example, Fig. || of paper Q 

2 )Ref. also the review article |jj and the papers listing in it 
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by the pseudogap decreases. The prediction discribed above is consistent with the 
experiment ||18|| 2 \ 



5. Concluding remarks 

It is of interest to note that the situation is entirely different if t' = 0. For in this 
case, Eq. ( |T2| ) reduces to 



surface, according to fllBD, is constant. 

Now, it is clearly that the peculiarity of the band structure of the cuprate plays an 
important role in understanding the pseudogap phenomenon in underdoped cuprate. 
This is the reason why our mean field SDW analysis of the pseudogap, based on the 
t — t' — U Hubbard model, meets with success. 

The mean-field solution has an antiferromagnetic long-range order. At sufficient 
doping concentration, the spin long-range order will be removed by fluctuations but 
there are still short-range orderings. We assume implicitly in our theory that the 
pseudogap structure, at least near the saddle point (ir, 0), is not sensitive to the 
long-range order and will survive in underdoped region, leading to the pseudogap 
phenomenon. 



One of the authors (H. S. Wu) would like to thank Prof. Z. Y. Weng for very 
valuable discussion. 




(13) 



It is in contradiction with the experiment P],rS], for the pseudogap along the Fermi 
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Figure Captions: 

Fig. [I]. The Brillouin zone and energy contour lines: The T point is at the middle 
of the Brillouin zone and the M points [(±7r, 0) and (0, ±7r)] are midway along 
the edges. The curves are the energy contour lines (t'/t = —0.16) with energy 
eik/i = —1.59, —0.64, —0.49 and —0.4, which are from inside to outside. 

Fig. H]. The period Brillouin zone: The solid curves are the energy contour lines with 
£k = 0. When the region II is shifted by the vector Q = (it, it), it coincides with 
the region I. In regions I + 1', > 0. In the region II, < 0. The region I' is 
called as the necklace region. 

Fig. |3|. The Fermi surfaces (t'/t = —0.16) in the quarter of Brillouin zone: The light 
curve represents the Fermi surface. 1 and 2 for the overdoping. 3 and 4 for the 
underdoping. The heavy dashed and the solid curve represent the necklace region 
boundary. 

Fig. [|. Our choice of the Brillouin zone. The heavy rectangle is our choice of the 
Brillouin zone boundary and the origin is at the M. The heavy square is the 
magnetic Brillouin zone boundary. 

Fig. |5|. This figure is the part of the magnetic Brillouin zone of the Fig. 4. The light 
curve represents the Fermi surface for the underdoping region (t'/t = —0.16). The 
dashed and the solid curve represent the necklace region boundary. The heavy 
solid lines are the magnetic Brillouin zone boundary. k x - and k y -axis are parallel 
with iWT and MX, respectively. <fi = aicta.n(k x /k y ). cf)' = a.rcta,n(k x / (it — k y )). 

Fig. |6]. The angle dependence of the pseudogap located at 0', Apg((f>'), for the hole 
doping concentration x=0.13 (t'/t = —0.18 and U/t = 0.8). The A<p' is the region 
where the values of Aps(4>') are all smaller than 2 meV. 
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Fig. |7|. The hole doping concentration dependence of the pseudogap located at 0' = 0, 
Aps(0), for t'/t = -0.18 and U/t = 0.8. The x indicates the hole doping 
concentration. The solid curve represents the pseudogap in the underdoping 
region. 

Fig. [§[ The hole doping concentration dependence of the region where the values of 
Aps(4>') are all smaller than 2 meV, A0', for the underdoping region (t'/t = —0.18 
and U/t = 0.8). The x indicates the hole doping concentration. 
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